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We study the first passage time processes of tiie anomalous diffusion on tlie self similar curves 
in two dimensions. The scaling properties of the mean square displacement and mean first passage 
time of the fractional Brownian motion and subordinated walk on the different fractal curves (loop 
erased random walk, harmonic explorer and percolation front) are derived. We also define natural 
parametrized subordinated Schramm Loewner evolution (NS-SLE) as a mathematical tool that can 
model diffusion on fractal curves. The scaling properties of the mean square displacement and mean 
first passage time for NS-SLE are obtained by numerical means. 



I. INTRODUCTION 

The anomalous or non-Fickian diffusive transports 
have attracted a lot of interest in the past few years. 
There is a wide range of heterogeneous or pre-asymptotic 
systems in the fields of physics, astronomy, biology, chem- 
istry, and economics, where anomalous diffusion occur 
This phenomenon is observed in diffusion on frac- 
tal structures in geophysical and geological media [I, Q , 
charge transport in disordered and amorphous semicon- 
ductors (4l, acceleration of particles inside a turbulent 
medium [5|, transport process in the biological systems 
0, and many other examples, for an extensive collection 
of references see [l| . The well-known examples of anoma- 
lous diffusions are the continuous time random walk (sub- 
diffusive processes) 0] and Levy flight [1| . In this type of 
processes, the mean square displacement (MSD) obeys a 
power law equation with respect to the time with ex- 
ponent < v < 2. Anomalous transport, especially 
the continuous time random walk and Levy flight can 
be studied within the fractional Fokker Planck equation 
approach Q. 

One of the most important characteristics in normal 
and anomalous diffusions is the first passage time (FPT) , 
which is defined as the time needed the dynamic variable 
to cross a given threshold value for the first time [lol . [Tl| . 
The FPT has been used to characterize diffusive pro- 
cesses in various systems such as the spreading of disease 
p^ . the passage of polymers and DNA in sub-diffusive 
media and membranes [l^jthe firing of neurons (l^ . 
animals searching for food [l5| and Levy stable random 
motion 

One of the interesting subjects of the diffusion problem 
is the study of FPT processes in fractal geometries such 
as percolating fronts, crack patterns, polymer chains, 
lightening paths, etc [l6l - [20j . In this work we are in- 
terested in simple fractal objects with fractal dimension 
I < df < 2 without any branch point. Since the coordi- 
nation number of all the points on the fractal is two it is 
easy to conclude that there should be lots of similarities 
between FPT processes in these systems and one dimen- 
sional systems. This was already discussed in [2l| and 
the important rule of the length of the fractal objects in 
FPT processes were understood. In [2l| the connection 



with the Schramm Loewner evolution (SLE) |23| were 
also discussed, especially the important rule of the defi- 
nition of the length in the SLE studies [U, were em- 
phasized. In this work we will generalize the work done 
in [2l|, in many different directions. To have an idea 
of fractal objects with fractal dimension 1 < df < 2 we 
will study the scaling exponents of MSD and mean FPT 
(MFPT) of the diffusing particles on loop-erased random 
walk (LERW), harmonic explorer (HE) and percolation 
fronts (PF) on the upper half plane. We will study differ- 
ent random walkers such as fractional Brownian motion 
and subordinated walk on the fractal curves. The scal- 
ing properties of the MSD of the walker and MFPT will 
be discussed. Finally we give a novel method to study 
diffusion on fractal curves by using SLE. We show that 
all the scaling behaviors discussed for the discrete fractal 
curves can be rederived by using subordinated natural 
SLE which is the time changed Schramm Loewner evo- 
lution. 

The central aim of this work is to use the subordinated 
natural SLE to find a new connection between diffusion 
on the self similar traces and SLE as a growth process. 
Our analysis of MSD and MFPT are the key points of 
this connection. 

The paper is organized as follows: In the next section we 
will fix the notation and introduce the scaling relations 
for the diffusion problem on scmi-lD fractal paths, where 
we measure the scaling exponents of MSD and MFPT, 
for two sided diffusion and diffusion with waiting times 
on the self similar curves with fractal dimension df. In 
the third section we use the Schramm Loewner evolution 
(SLE) as a mathematical model to introduce new classes 
of diffusion processes (subordinated SLE). The results of 
this section are compatible with the two-sided diffusion 
and diffusion with the waiting time on the discrete fractal 
paths. In the last section we conclude our findings. To 
be self explanatory we add three appendices explaining 
the details of our simulation methods. 
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II. FIRST PASSAGE TIME IN LATTICE 
FRACTAL INTERFACES 

We begin by considering the diffusion problem on semi- 
ID random curves. For our purposes we restrict our- 
selves to the fractal curves with HausdorjJ dimensions 
1 < d/ < 2 [l^ that start from the origin and remain 
in the upper half-plane. To understand statistical prop- 
erties of diffusion along such inhomogeneous paths it is 
important to first introduce diffusion problem in ID case. 

A. First passage time statistics in one dimensional 
domain 

Consider a one-dimensional diffusion Xt with dynamics 
dXt = a{Xt)dt + adWt, (1) 

where Wt is a one dimensional stochastic process. The 
interval for the solution of Eq. ([TJ is defined as closed on 
the left hand side Xa = and open on the right hand side 
Xb = R- These special choices force diffusing particle to 
move only in a; > 0. 

An interesting problem in the theory of stochastic pro- 
cesses is finding the time that a particle reaches a certain 
level. The problem of finding this time is called first pas- 
sage time (FPT) (lol . IT]| . The first passage time is the 
time taken for the process having started from x — Q 
to be reached in x ~ r [lij 

Tr = inf{t > 0\Xt = r}, (2) 

where the infimum for every subset S of real numbers is 
denoted by inf{S'} and is defined to be the biggest real 
number that is smaller than or equal to every number in 
S. 

Clearly FPT is a random variable which varies from 
one sample of Xt to another one. In general we are in- 
terested in those processes where the two statistical vari- 
ables {Xt) and {Tr) have scaling behavior 

OC t\ {Tr) cx r^. (3) 

One example of Eq. ([T])isa = 0,(T = l and Wt — 
\B^ I, where B^^ is the fractional Brownian motion (fBm) 
process. Note that, fBm with Hurst index < iJ < 1 is 
the only self-similar Gaussian process with stationary of 
increments [2^ . The correlation function of fBm is 

{B"{t)B"{s))^[\tr + \s\^'' -\t-s\'"]. (4) 

Absolute value of B^ is used to force diffusing particle 
to move only in the region a; > 0. The particle's position 
Xt = \Bf\ is a stochastic variable with 

{X^) = {\B^t) ^t^". (5) 

The scaling relation in Eq. ([S]) is in agreement with Eq. 
^ and shows that v = 2H [23] . Diffusion is said to 



be anomalous if 7^ 1, where 0<i/<l(l<i/<2) 
indicates the sub-diffusive (super-diffusive) behavior. 

It is easy to determine the exponent (3 = 1/H ana- 
lytically for the mentioned boundary condition. To find 
the scaling parameter (3, consider a random process Xt 
as Xt/x^/H = j^t- Now one can replace r — > Ar in Eq. 
([2]) to find Txr- It is straightforward to show that Txr is 
the same as X^/^Tr in the distributional sense. By this 
scaling argument, one can observe that MFPT is given 
by 

{Tr) OC r^/". (6) 

It is worth mentioning that by choosing Wt = \Bt\, where 
Bt is a Brownian motion (fBm with H = one can use 
the Fokker-Planck (FP) equation, which describes the 
space-time evolution of the probability density function 
(pdf) of Xt, to find FPT distribution. More details can 
be found in 0. 

Although studying diffusion in one dimensional sys- 
tems is interesting for its own sake there are many ex- 
amples of diffusion in self-similar interfaces with fractal 
dimensions 1 < dj < 2. We will generalize the above 
arguments to self-similar interfaces with arbitrary fractal 
dimension and will study the statistical properties of the 
diffusion on the fractal curves. 

To motivate our method of measurement of the scal- 
ing parameters v, f3 for diffusing particle on the fractal 
curves we first consider models on a lattice domain, e.g., 
the loop-erased random walk, harmonic explorer and the 
percolation explorer process. 

a. The loop- erased random walk 

The LERW on the square lattice domain is a random 
walk with erased loops when they appear. This process 
is stopped when it reaches a given length. To produce 
LERW curves started from the origin and conditioned to 
be in the upper half plane one can use reflecting boundary 
condition on the real axis for the random walker. The 
fractal dimension of LERW is 5/4 [131 . 

b. An explorer processes 

Explorer processes (EP) on honeycomb lattice such as 
percolation front with df = ^ and harmonic explorer 
with df = I are used as other classes of fractal interfaces 
[2^ [25| . To construct an EP path with a fixed number 
of steps N, we used from a class of explorer processes 
on the honeycomb lattice. This process is named over- 
ruled harmonic explorer [sij . For numerical analysis, we 
simulated this process on the extremely large rectangular 
domain, where it can approximate the upper half plane 
(see appendix. A). 

Although other studies such as HJ] have also used 
self similar traces to study the FPT problem they have 
not characterized the scaling relations in the measurable 
quantities such as MSD and MFPT for the general diffu- 
sion processes e.g. two sided diffusion and diffusion with 
the waiting time which we have studied in our simula- 
tions. 
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B. Two sided Diffusion on the fractal paths 

In this subsection wc study the statistical properties 
of the diffusing particles along self-similar curves. An 
interesting problem in this direction is the determination 
of the scaling exponents of the random displacement. 

For one dimensional domain with reflecting boundary 
condition, we mentioned in section 2.1 that the MSD and 
the MFPT according to Eq. ([5]) and Eq. (|6]), obey scaling 
laws with the exponents v = 27J and (3 ^ 1/H. In the 
following we explain how to use the discrete random walk 
model to simulate stochastic process Xt on ID domain 
with reflecting boundary condition on a; = and we then 
favourably apply this model to random process on the 
discrete fractal curves. 

First we consider the random walker on the one dimen- 
sional discrete lattice. This random walker started from 
the position x = Q ai = Q and at the time f„ = n5t 
moves one step to right (left) when {l-B^^J — \B^\} > 
i{\B"+st\ - l-Sf 1} < 0). The normal random walk 
(discrete version of Brownian motion) corresponds to 
H = 1/2. This random process corresponds to two sided 
diffusion on ID domain. Following the idea presented in 
one dimension, we obtain the statistics of two sided dif- 
fusion on the fractal curves. To this aim we consider a 
random walker with position coordinates Xn and Yn for 
the n*'* walk where {Xq, Yq) = (0, 0) is the start position 
and it moves back and forth along the discrete self-similar 
curve. In order to simulate random walker on the curve 
started from the origin and remains in the upper half 
plane we used from fBm process \B^\. Using this cor- 
related stochastic process wc define another stochastic 
process so that S'^^^i ^ + 1 (S'^_^i = - 1) 
when - |Sf 1} > ({|i?,%| - |i?f |} < 0) with 

the initial value Sq = 0. The random walk position can 
be defined by Xn = x{Sn) and Yn = y{Sn) where x{i) 
and y{i) are the position components of the i*'' point of 
the curve. 

We now study numerically scaling dependence of 
(i?^) = {Xn + Yn) and (r^) to n and r for many random 
walkers moving along such self similar one-dimensional 
objects. Especially we study the scaling forms of {R^-^ 
and (r^). 

To study the scaling laws in two sided diffusion, it 
should be mentioned here that the scaling exponents 
v{H, df) and /3(i?, d/) in two sided diffusion on the frac- 
tal curves are in general a function of Hurst parameter H 
and geometrical dimension df. Within the EFT statis- 
tics of one dimensional diffusing particle approach, one 
can consider the self-similar curve as a one dimensional 
non-straight line with length /. For this semi-ID object 
using Eq. ([6]) one can observe that (t;) oc V-^^ . On the 
other hand for a fractal curve (with length I) inside a 
circle there is a scaling law / oc r'^^ where r is the radius 
of the circle. Under these assumptions, the MFPT reads 

as {Tr) (X r"^!/" . 

The very same method can be applied to scaling law 



of MSD. Same as before we consider the fractal curve as 
a semi-ID object. There is a scaling relation (/^) (x 
(similar to Eq. [S|) for the position of diffusing particle 
In after n walks along such semi-ID curve. In addition, 
the scaling relation Z„ ~ ufj for the fractal curve is well- 
known, where i?„ is the radius of the semi-circle enclosed 
Z„-th walks. We therefore obtain a universal law c>c 
^2H/df Jqj. j^gj3 gjjgj diffusin g particles on the 

fractal curves. Using these arguments we expect 

v{H, df) - 2H/df, p{H, df) = df/H, (7) 

These equations are in agreement with the scaling ex- 
ponents in the Eqs. ([5]) and (O in the c?/ — >■ 1 limit. 
The above results for H = ^ recover the predictions in 
[2H . Table |T] summarizes our numerical results for the 
two scaling exponents v and /?, where it shows that our 
results are well compatible with the predictions in Eq. 
([7]) for two sided diffusion on the fractal curves. In our 
measurements we used from 50000 fractal curves and 10 
independent realizations of random process Sn per curve 
for each numerical test. 



Model iy{H ^ 0.8) I3{H = 0.8) u{H = 0.9) /3(iJ = 0.9) 

LERW 1.280 ± 0.001 1.56 ± 0.01 1.437 ± 0.005 1.385 ± 0.005 

HE 1.066 ±0.001 1.88 ±0.01 1.195 ±0.005 1.670 ±0.003 

PF 0.92 ±0.01 2.17 ±0.01 1.029 ±0.003 1.945 ±0.005 



TABLE I: Numerical values of the scaling exponents v and 
l3 for the two sided diffusion on the fractal curves. To mea- 
sure these exponents we used from correlated process with 
H=0.8 and 0.9. 



C. Diffusion with waiting times on the fractal 
curves 

In this subsection we study anomalous motion of a free 
particle with waiting time on the self-similar curve. In 
the preceding section, we have obtained statistics of walk- 
ers moving back and forth randomly along self-similar 
discrete curves. The physical time between two consecu- 
tive steps of walks in the previous examples is equal to a 
constant At. 

We can also consider the random time elapsing be- 
tween two consecutive jumps of a diffusing particle. The 
particle starts from the origin and trapped in site n for 
some random time. These positive random waiting times 
T„ are identically distributed random variables each hav- 
ing the same probability density function ^'(t) 0,01 • The 
role of the waiting time forces us to identify operational 
time S't" = J2n '^n and physical time t„ = nAt. The 
physical time n is always accelerated against the strictly 
increasing random time S^^ . This is so, the random time 
S't" is called subordinator. As mentioned in (32| - [3^ it is 
described as 
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S? 



^f{T„l ■■ U{Tm) > tn}, 



(8) 



where U{Tm) is a-stable subordinator (0 < a < 1) and 
Tm = toAt. The random process Sf^ is cahed the inverse- 
time Qf-stable subordinator. The above process has nei- 
ther stationary nor independent increments but it is easy 
to show that we have distributional scahng 5""^ = X^Sf 
which leads us to the following symmetry for the subor- 
dinated Brownian motion 



B{sr). 



(9) 



Although B{S^) is self similar with Hurst exponent a/2 
it is not fractional Brownian motion because it docs not 
have a Gaussian distribution and it does not have sta- 
tionary increments. The process S"" is strictly increasing 
and it tends to infinity for t — > oo and so it is a good pro- 
cess to consider as the time, whereas in the a — >■ 1 limit 
the subordinated time converges to the physical time. In 
our study we generated the subordinator S*"^ following 
the reference [33] (see appendix. B). 

We study statistical properties of position coordinates 
Xn = x{Sf ) and Yn = y{Sf ) for a subordinated walker 
with operational time Sf^ moves on the discrete self- 
similar curve with fractal dimension dj and position com- 
ponents x{i) and y{i). For this class of particle diffusion 
in fractal path, we again expect universal scaling depen- 
dence of the MSD [{Rl) = {X^ + Y^)) and the MFPT 
{{tj.)) to the geometrical parameters n and r as 



v{a, df) = 2a/df, 



I3{a, df) ~ df /a. 



(10) 



Our numerical results for a = 0.8 and a = 0.9 (see Table 
HI)) arc in good agreement with the scaling exponents in 

Eq. (nni). 



used three statistical models, loop-erased random walk, 
harmonic explorer and percolation interfaces on the lat- 
tice. 

The scaling limit of the lattice models as the lattice 
spacing goes to zero corresponds to Schramm Loewner 
evolution. This mathematical model is defined in the 
complex plane and it was introduced by Schramm [23| . 
SLE is based on the Loewner equation 



dtgt{z) 



(11) 



where the real-valued function is called driving (or 
forcing) function, which determines all the properties of 
SLE. Loewner showd that for any non-intersecting curve 
parametrized by a complex function "f{t) in the upper- 
half plane H. there exist a conformal map gt{z), which 
maps upper half plane minus curve and the region which, 
is separated from infinity by the curve {hull: Kt) W\Kt^ 
to the upper half plane M |35|. 

Ordinary SLE is the Loewner evolution with ^4 = 
^/nBt , where Bt is the Brownian motion with mean zero 
and E\BtBs] = min{t, s) and also with diffusion con- 
stant K > [23j. These properties ensure that the curve 



is conformally invariant. 

SLEk is a random conformally invariant curve with the 
fractal dimension d/ = 1 + k/8 (0 < k < 8) [13, HI. The 
scaling limits of LERW, HE and PF are SLE^ with k = 
2, 4 and 6, respectively [lillO,!!^- 

The SLE„ curve (7(i)) is parametrized with the time 
t. On the other hand the lattice models (LERW, HE 
and PF) usually have a natural parametrization given by 
the number of steps with equal length along the curve. 
In general the scaling limit of the lattice models are 
not the same as SLE^, where this difference comes from 
parametrization of each model (38| . We will use from an 
appropriate method to re-parametrize the SLE^ curve. 



Model v{a = 0.8) /3(q = 0.8) v{a = 0.9) I3{a = 0.9) 

LERW L280 ± 0.001 L56 ± 0.01 1.440 ±0.001 1.39 ±0.01 

HE 1.065 ±0.005 1.88 ±0.01 1.200 ±0.001 1.67 ±0.01 

PF 0.920 ±0.005 2.18 ±0.01 1.030 ± 0.005 1.95 ±0.01 



TABLE II: Numerical values of the scaling exponents v and /3 
for subordinated diffusion on the fractal curves. To measure 
these exponents we used from subordination with a=0.8 
and 0.9. 



III. FIRST PASSAGE TIME AND SCHRAMM 
LOEWNER EVOLUTION 

In the preceding section we studied scaling exponents 
of MSD and MFPT for some important examples of two 
sided and subordinated diffusion on the discrete self- 
similar curves in the upper half plane. In our study we 



A. Natural parametrized SLE 

First step in simulating SLE is time-step discretiza- 
tion. Let us consider a partition of the time interval 
[0, t], where it is discretized into = to < ti < t2 < ■ ■ ■ < 
tn = t. One method to simulate SLE is the foregoing ap- 
proximation with the equally spaced discrete time points 
ti = idt. In this method the points Zi on the curve 'y{t) 
is given by an iteration process Zi = /i o /2 o • • • o (^^ ) , 

where fj{z) — \J [z — ^j)'^ — 4dt + is the inverse con- 
formal map and is the discretized drift, where it will 
be approximated by a piece- wise constant function in the 
uniform partition of the time interval [{i — l)dt, idt]. No- 
tice that the conformal map fi{^i) can produce a small 
slit at £,j with length Li ~ lm(fi{£^i)) = 2y'ti+i — ti on 
the upper-half plane. In this method the two-dimensional 
distances Z,; =| 7(<i)— 7(ti_i) | are extremely non-uniform 

We hereby, require the natural parametrized SLE^ (N- 
SLE^) curve, where it is the discrete SLE^ curve {7^} 
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with an approximately equal step length | 7(ti)— 7(ti_i) | 
~ A. There are some mathematical and numerical pro- 
cedures used to find a sensible definition of N-SLE„ 
[UliSlill (see appendix C). 



B. Subordinated SLE 

In order to understand the scaling relations for the 
models of diffusion with waiting time on the self-similar 
traces and also two sided diffusion on the fractal curves, 
we present here the subordinated version of SLE^ . 

The motivation of our approach comes from this idea 
that the probability distribution of the point at the tip 
of the SLE^ trace satisfies Fokker-Planck equation (FPE) 
[39} which basically one can think about it as the FPE 
of the position of a particle in the fractal interface. In a 
similar way one may think about the FPE for the prob- 
ability distribution of the tip of the N-SLE and also sub- 
ordinated N-SLE (NS-SLE) curve in the continuum limit 
and study diffusion. In principle it should be possible to 
calculate analytically MSD and MFPT for these semi-ID 
interfaces by using the FPE in two dimensions. Unfor- 
tunately we do not know how to write the FPE of the 
N-SLE and NS-SLE, therfore we just calculate numeri- 
cally the scaling properties of the tip of NS-SLE as the 
diffusion process on the fractal paths. We will show that 
the scaling behaviors of the subordinated N-SLE are sim- 
ilar to the lattice models. 

For normal SLE given by Eq. pT|) the time variable 
is deterministic but we would like to set this variable 
as an internal parameter t which is also stochastic and 
strictly non-decreasing, this is called subordinating the 
process by the inverse time a-stable subordinator (see 
appendix. B). 

Using the above definition one can consider Loewner's 
map with the new time as gsp (z) which is still scale 
invariant in the following sense: the conformal map 
~9s?{z) = ^5sj,(A"/2z) with B{Sr) := jh^BiS^,) 
satisfies the same Loewner equation as gt{z). The above 
scale invariance enforces scale invariance of the curve. 

We have discussed the simulation of SLE and N-SLE in 
the last section. The simulation of subordinated SLE (S- 
SLE) and natural parametrized subordinated SLE (NS- 
SLE) are similar. The only difference is the confor- 

mal map f,{z) ^ ^ [z ~ ^S^^))^ ~ AdS'^^ + ^S?,) where 

dS^^ = Sf^ - Sf^_^ and ^ to < ti < t2 < ■ ■ ■ < tn = t. 

The time steps Aj = dt {ti = t^^*^ *^^sc 

of S-SLE^ are selected uniformly and in the NS-SLE^ 
the non-uniform time steps Ai are computed by using 
Jacobian scheme. The only difference between S-SLE" 
(NS-SLE^) and normal SLE«; (N-SLE„) is in the growth 
process of each of them. In the first case the tip of the 
curve has waiting time according to the a-stablc Levy 
process. 

As discussed earlier, the scaling exponents v and /? for 
subordinated walk along discrete fractal interfaces arc 



defined explicitly in Eq. (jlO[) where they are in agree- 
ment with numerical simulations. We will consider a tip 
of NS-SLE" as a subordinated growth process where it is 
a mathematical model for subordinated walk along frac- 
tal curves. The scaling exponents of this subordinated 
process is collected in Table. Illli where they are in a 
good agreement with Eq. (|10p and also with the numer- 
ical simulation of the subordinated random walk along 
the self-similar discrete curves (see Table. HI)) . 



NS-SLEg v{a^Q.9,) /3(q = 0.8) ;/(q = 0.9) ^(q = 0.9) 

K = 2.0 1.28 ± 0.02 L52 ± 0.05 L42 ± 0.04 L35 ± 0.05 

K = 4.0 1.02 ± 0.05 1.84 ± 0.04 1.16 ± 0.05 1.63 ± 0.05 

K = 6.0 0.92 ± 0.03 2.15 ± 0.05 1.00 ± 0.03 1.90 ± 0.06 



TABLE III: Numerical values of the scaling exponents v and 
P for the tip of the NS-SLE° curves with k =2.0, 4.0 and 6.0 
and a=0.8 and 0.9. 

We also notice that another way to subordinate the 
forcing function in Eq. pT|) is based on the iterated 
Brownian motion (42j . Consider two stochastic pro- 
cesses Bt and , where the first one is the Brownian 
motion and the second one is the fractional Brownian 
process. The iterated Brownian process is defined as 
B{\Y^\), where \Y^\ corresponds to the non- negative 
random time. It is easy to verify that the fractional 
Brownian time Brownian motion _B(|yj^ |) is a self-similar 
process of index H /2, that is, for any A 

i?(|y,f I) = B{X"\Y,"\) = X"/^B{\Y,"\). (12) 

Simulation of the the natural parametrized version 
of fractional Brownian time SLE (NF-SLE^) IS simi- 
lar to the NS-SLE case. First consider discrete times 
Q = to < ti < t2 < ■ ■ ■ < tn = t and the confor- 

mal map f,{z) = y^(z - ^(|r,f |))2 - 4d|y,f | -f ^(|y,f |), 
where the infinitesimal values of the local time cZlFj^l = 
\Y^\ — \Y^__^ I can get positive and negative values. The 
length of NF-SLE^ curves increases for dl^t^l > and 
decreases for c?|yi^| < 0. This dynamical process is very 
similar to the two sided diffusion on the lattice fractal 
models. Our estimations for the two scaling parameters 
V and (3 for MSD and MFPT, Table [iVl are in a good 
agreement with the predicted values in Eq. ^ and also 
numerical results coming from the lattice models (see Ta- 
ble [J). 



IV. CONCLUSION 

To conclude, we studied the diffusive dynamics of the 
random processes on the self-similar curves and measured 
the scaling exponents of mean squared displacement and 
mean first passage time expressed in Eq. ([3]). The various 
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NF-SLE^ iy(H ^ 0.8) P{H = 0.8) u{H = 0.9) P{H = 0.9) 



K 


= 2.0 


1.28 ±0.01 


1.57 ±0.01 


1.44 ±0.01 


1.40 ±0.02 


K 


= 4.0 


1.05 ± 0.03 


1.90 ± 0.02 


1.16 ±0.05 


1.70 ±0.04 


K 


= 6.0 


0.90 ± 0.03 


2.20 ±0.02 


1.01 ±0.03 


2.00 ±0.05 



TABLE IV: Numerical values of the scaling exponents i/ and 
P for the tip of the NF-SLE^ curves with k =2.0, 4.0 and 6.0 
and H=0.8 and 0.9. 

scaling exponents for MSD and MFPT are obtained nu- 
merically for two sided diffusion and diffusion with wait- 
ing time on three discrete fractal curves, i.e. loop-erased 
random walk, harmonic explorer and percolation front. 
It appears that the exponents only depend on the fractal 
dimension d/ of the curves and the scaling exponent H 
for the two sided diffusion and a for the subordinated 
diffusion. 

Finally, we rederived the properties of the anomalous 
diffusion (FPT, MSD) on the discrete fractal curves with 
subordinated version of the natural parametrized SLE. 
Our results offer a new method to investigate diffusion 
in the fractal interfaces. We believe that these results is 
a starting point for the development of the subordinated 
version of SLE. 



Appendix. A 

To find harmonic explorer and percolation front, we 
used overruled harmonic explorer process on the very 
large rectangular domain. This domain on the upper half 
plane as shown in Fig. 4 is splitting into three parts, a 
left boundary with yellow condition and a right boundary 
with blue and also uncoloured inner part. This bound- 
ary condition is used to limit the EP path to this part of 
half-plane to start from rg and stop when reach to or 
length N. The explorer process is the unique path from 
the origin. In each step, there is a yellow hexagon on the 
left and blue one on the right [25| . 

To generate this path dynamically, a growth process 
starts from the point tq on the lower boundary. In the 
first step the colour of face /i in front of tq is chosen so 
that to make it blue or yellow and the explorer is forced 
to turn left or right, respectively. To choose colour of face 
/i, a random walker starts from /i and it stops when it 
crosses the rectangle's boundary for the first time. Now, 
the colour of /i with probability < p < 1 is yellow if 
the touch boundary is yellow. Note that two stochas- 
tic operations are used to colour one hexagon, a random 
walker to find boundary colour and a coin to does or does 
not accept the boundary colour. The new tip of explorer 
path is located in the position ri and a new face (/2) 
should be coloured with the same restriction. In par- 
ticular, the outcome of explorer process with p = i,l 
as shown in Fig. 5 are percolation front and harmonic 
explorer, respectively. The fractal dimension of the over- 




FIG. 1: A rectangular domain with appropriate boundary 
conditions is used to build an explorer path. Top: Three parts 
of the domain which consist of left (right) boundary with 
yellow (blue) hexagons and uncoloured hexagons. Middle: 
First step to identify colour of face /i. A random walker 
moves around uncoloured sites to hit the boundary. In this 
example, the colour of face fi with probability p (1 — p) will 
be blue (yellow). Notice that walker turn to right (left) when 
the yellow (blue) is selected. Bottom: A complete exploration 
process in a rectangle. 



ruled harmonic explorers has a linear relationship with p 
and it is conjectured to be d/ = 2 — | [3l| . 

Appendix. B 

The celebrated subordinated random time S'f , is given 
by the Eq. ^ where it can be efficiently generated by the 
algorithm proposed in (33| . As mentioned earlier the a- 
stable subordinator U(t) denotes the strictly increasing 

Levy motion with Laplace transform (e^'^^''^-') = e^'^*^ 
[s^ ]. The first step in simulating S" begins with approx- 
imating of the strictly increasing a-stable Levy motion 
U{t) on the discrete times t; = iAr {i ~ 0,1,..., A/). 
The numerical integration of the process U{t) for < 
a < 1 yields 

C/(r,+i) = U{rj) + At1/"L„(/3), (13) 
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K = 2.0 



K = 4.0 



K = 6.0 




FIG. 2: The overruled harmonic explorer path with length 
A'' = 3 X lO"*, Left: Harmonic explorer with p = 1.0 and 
fractal dimension dj = 3/2, Right: Percolation explorer with 
p = 1/2 and fractal dimension d/ = 7/4. 



where La{P) is a Levy stable random variable with pa- 
rameter /3 and U{0) = 0. We use skewed Levy-stable 
distribution (/3 = 1), to ensure U{t) gets almost increas- 
ing random process [34| . It can be generated by 



L„(1.0) 



sin[«(y + f )] ^ cos[V - a{V + |)] 



[cos(\/)]i/" 



x{- 



ly 



(14) 

where y is a random variable with uniform distribution 
between (— -f, f) and W has exponential distribution 
with mean 1. For the time horizon T, the summation 
process in Eq. ends when we get U{tm-i) < T < 

U(tm)- One can observe that U{t) is strictly increasing 
and M always exists [11] . 

Now, for every t,; € (0 = io < < ^2 • • • < = T), 
we find Tj such that C^(tj_i) < ti < U{Tj), and from the 
definition in the Eq. ([5]), we can define S" = tj. From 
Eqs. (HH) and ^ it is clear that ^^(l) = i and S^. = U 
in the a — > 1 limit, where at this limit subordinated time 
converges to the normal time. 



FIG. 3: The N-SLE« curves with length iV = 1 x 10*, A = 
0.001 and k=2.0, 4.0 and 6.0 from left to right. 



Note that computing points along the N-SLE^ curve 
requires this adaptive choice of and the total time 
with these non-uniform time steps will be equal to = 
Y^l^^i Ai. In this procedure the distances between two 
sequential points U k. \ approximately remain constant. 

In our study we followed one straightforward moti- 
vation for computing the Jacobian. If one consider 
hi{z) = fi{z + — The conformal map /ii(0), 

maps the upper half plane onto the upper half plane plus 
a slit. The length of this slit equals to 2y/Ki and the po- 
sition on the real line equals to 6i = S,i — Ci-i • Following 
a simpler strategy one can decompose the incremental 
map hi{z) to h.i{z) = Ts^ o where = Vz'^ - 4Aj is 
the slit map and Tg. (z) = z + Si is a translation map by 
the real value Si. The z*'* points of the the SLE or N-SLE 
curve computed from j{ti) = 5„(0) where 

5„(z) = o 0f o T^, o o . . . o Ta. o 0f (z). (17) 

We now consider a new format of Eq. as 

A 



I gliio) 



(18) 



Appendix. C 

A standard procedure [i^l to find the half plane N- 
SLE^ trace is based on a change in the size of the z*'' slit 
length Li (It is a function of time step parameter A.^ = 
ti — ti-i as Li = 2\/Ai), by the Jacobian. The Jacobian 
I Ji-i l~l - ^i-i)G-'_i(^i_i) I of the conformal map 
G; = /i o /2 o • • • o /i acts on the corresponding segment 
to rescale the length Li for the i*'' slit by 



L,. 



A 



J 71- 



(15) 



where A > is the step length. For a piece-wise constant 
Brownian process = ± \/ nAi (the sign of a/kAJ" 
is chosen randomly according to the uniform probability 
distribution), the above approximation yields 



A 



2V^iG:'_i(e.-i)r 



(16) 



where 

n-2 

^'(0) =1 0;;(o) I n i'/'n-i-,(r,) 

In the above equation a F,- is defined as 



(19) 



r, = r5„_,o</,H_^.oT,_^.^,o0H_^.+iO. . •or,„o</,H(0). (20) 

Following [40| the proposed method of approximating 
sample paths of N-SLE^ {^{ti),i = 0,...,N}, consists 
of six steps: (1) Set the constants A, k and N. (2) Set 
n = 1 and Ai = 1. (3) Compute V kA„ according to 
steps 1 and 2 with a random sign (±) with equal prob- 
ability. (4) Calculate 7(tn) = 5n(0) using the iteration 
map as we said in Eq. p?)) . (5) Compute A„+i using 
Eq. (HH), (HH) and (6) If n < iV increase n by 

one and repeat steps 3 to 6. The typical curves of the 
N-SLE„ for K =2, 4 and 6 are presented in Fig. 6. 
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